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Abstract 

Some effects of irradiation on time-dependent accretion-disk models for black hole X-ray novae are 
presented. Two types of irradiation are considered: direct irradiation from the inner hot disk and indirect 
irradiation as might be reflected by a corona or chromosphere above the disk. The shadowing effect of the 
time-dependent evolution of the disk height and consequent blocking of the outer disk by the inner and 
middle portions of the disk from the direct irradiation is included. The direct irradiation of the disk by 
inner layers where the soft X-ray flux is generated is found to have only a small effect on the outer disk 
because of shadowing. Mild indirect irradiation that flattens, but otherwise does not affect the light curve 
substantially, still has interesting non-linear effects on the structure of the disk as heating and cooling waves 
propagate. The irradiated disks do not always make simple transitions between the hot and cold states, 
but can linger at intermediate temperatures or even return temporarily to the hot state, depending on the 
irradiation and the activity in adjacent portions of the disk. 
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1. Introduction 

There is now broad recognition that the most common 
form of black-hole binary candidates in the Galaxy are 
transients (Chen et al. 1997). The cause of the primary 
rise in outburst is most likely to be due to the disk ther- 
mal instability associated with the ionization of hydrogen 
and helium (Cannizzo et al. 1982; Mineshige, Wheeler 
1989). The disk-instability model can also produce ultra- 
soft X-rays (^ keV) and a strong correlation between the 
optical and X-rays in the outburst evolution (Mineshige 
et al. 1990a). There remain many interesting issues, 
however. The original disk- instability model (without ir- 
radiation) failed to produce a plateau in the ^-magnitude 
like that 100-200 days after the peak in some systems. 
Although the disk-instability model can to account for 
the beautiful exponential decays (Cannizzo et al. 1995), 
it has difficulty in producing the refiares (or secondary 
maxima) commonly observed in X-ray novae at 50-70 
days after the main peaks. These were suggested to be 
due to irradiation of the outer portions of the disk (Mi- 



neshige, Wheeler 1989; Mineshige 1994). 

We examine here some of the aspects of the effect of 
irradiation on the propagation of heating and cooling 
waves in time-dependent disks that are unstable to the 
ionization thermal instability. Preliminary versions of 
this work were presented by Kim et al. (1994), Wheeler 
et al. (1996), Kim et al. (1996), and Wheeler (1997, 
1999). A precursor to this work was that of Saito (1989), 
who was the first to consider time-dependent, irradiated 
disks with some simplifications (discussed later) . Because 
this is the first study of time-dependent irradiation in 
black-hole disks, we concentrate on qualitative effects in 
systems like A0620-00. The model assumptions, includ- 
ing a prescription for treating irradiation, are given in 
section 2. The global effects of irradiation are presented 
in section 3 and the effects of irradiation on the thermo- 
dynamics of the heating and cooling instability are given 
in section 4. A discussion and conclusions are given in 
section 5. 
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2. Disk-instability Model for Irradiated Accre- 
tion Disks 

2.1. Irradiation 

Two types of irradiation are considered: direct irradi- 
ation from the innermost hot disk and irradiation that 
might be reflected by a corona or disk atmosphere or 
chromosphere above the disk (Meyer, Meyer-Hofmeister 
1994). The X-ray luminosity of the irradiation is given 

by 



Lx{t) = eMinC^ 



(1) 



where the efficiency e is 0.057 and Mjn is the mass- 
accretion rate at the inner edge of the disk R^^, taken 
to be 3rg as for a Schwarzschild black hole. We adopt 
a simple model for indirect irradiation by assuming that 
the luminosity of equation (|^) effectively arises from a 
point at the center of the disk {R — 0). The indirect 
irradiation flux is then given by 



Fi{R,t)^C, 



Lxjt) 



(2) 



where Cx is a constant. More sophisticated models in- 
corporating radiative transfer would give Cx = CxiR^t) 
(see discussion in Tuchman et al. 1990), but such time- 
dependent radiative transfer is too difficult at this time. 

For direct irradiation, we follow the prescription of 
Fukue (1992; see also Kusaka et al. 1970). This model as- 
sumes that the irradiation from a hot, geometrically thin 
annulus near the inside of the disk can be approximated 
by that from an infinitesimally thin, filled, uniformly ra- 
diating surface centered on the black hole. For regions 
in the disk at distances much larger than the radius of 
the annulus, this leads to the following expression for the 
total flux incident on one surface of the disk: 



Fd{R,t) = il-A) 



Lxjt) d_ 
47rR ' dR 



H{t) 



R 



where A is the X-ray albedo. 

In practice, this irradiation is computed to arise from 
a disk with the radius of the inner zone and from a height 
above the disk plane equal to the disk height, iJ, of the 
innermost computed zone. The latter is a function of 
time and is used to compute the rays for direct irradiation 
of the outer disk. 

Equation (^) is equivalent to that adopted by King et 
al. (1997, see their equation 6). It is smaller by a fac- 
tor of H/R than the prescription used by van Paradijs 
(1996) for neutron-star systems. For simplicity, we take 
the X-ray albedo to be A^ 0.5. De Jong et al. (1996) 
have deduced an albedo of ^ 90 % for low-mass X-ray 
binaries, but the specific choice of the albedo does not ap- 
preciably affect the disk structure, since shadowing sup- 
presses direct irradiation (see subsection 3.1). Equation 



(^) shows that direct irradiation is a function of the gra- 
dient of the disk height, and is hence very sensitive to 
small variations in the disk profile. Analytic models as- 
suming d\nH/dlnR = constant may exaggerate the ir- 
radiation by minimizing the shadowing that may occur 
even in the steady state. They are certainly not adequate 
for time-dependent disks. The height of the region from 
which the irradiation arises varies with time and the disk 
height profile varies with both the irradiation and time. 
This makes the shadowing of the outer disk a complex, 
time-dependent phenomenon. We assume that no direct 
irradiation flux is added to the shadowed portions of the 
disk, while the indirect irradiation heats the whole disk 
in accord with equation (j^. 

2.2. Basic Equations and Boundary Conditions 

Once Mill, and hence Fi{R,t) and Fd{R,t), are deter- 
mined, the irradiated flux is added to the flux generated 
internally in the disk by viscous heating through an im- 
plicit numerical method, and other physical variables are 
computed in each time step (see Tuchman et al. 1990 
for details). A modification of the outer boundary con- 
dition by irradiation affects the whole vertical structure 
of the disk, including the mid-plane values (Tuchman et. 
al. 1990). This is a reasonable first approximation, but 
could be improved by including more sophisticated radia- 
tive transfer effects. In this formulation, the irradiation 
fluxes are parameterized by two constants, Cx and A, 
which are, by assumption, independent of the radius and 
time. 

The basic equations of the time-dependent thermally 
unstable accretion-disk model (for details, see Kim et al. 
1992 and references therein) are the mass conservation, 



(9S BM 
2nR— ^— + 2TrRSoiR,t), 



(4) 



(3) the angular momentum conservation; 



27ri?^ = ^iM^+2nRSo{R,t)iK-4^{2nR^WU5) 
at oR oR 

and the vertically averaged energy equation; 

¥ (I +^^4) = 0+s+^^d+Fi-0-,+Qdif,(6) 

where we use cylindrical coordinates {R,(j),z), with the 
origin at the black hole and the z-axis perpendicular to 
the plane of the disk {z = 0). Also, E is the surface 
density, So{R,t) is the mass source function defined by 



Mrit) = / 27rRSo{R,t)dR, 



(7) 



with Mt being the mass transfer rate from the compan- 
ion star, i = {CMrY^"^ is the specific angular momentum 
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of matter in the disk, £k is the specific angular momen- 
tum of the stream from a companion star, M is the mass 
of black hole, W {— j aPdz) is the integrated viscous 
stress (with a and P being the viscosity parameter and 
pressure, respectively), Tc is the central temperature of 
the disk, vr = — M/27ri?E is the radial velocity of the 
matter, Cp is the specific heat with constant pressure, 
and f2K is the Keplerian angular frequency. 

There are several terms on the right-hand side of equa- 
tion (^). The disk heating in the non-irradiated disk is 
given by the viscous heating; 







-wn. 



(8) 



The external-heating terms due to indirect and direct 
irradiation of the disk, and i^d, are given by equations 
(^ and (^). In the present time-dependent study, we 
included irradiation as heating in the vertically averaged 
structure, and did not solve the vertical structure of the 
disk with irradiation input from above the surface (cf. 
Tuchman et al. 1990). The radiative cooling is given by 



<3rad — '^^cffi 



(9) 



where a is the Stefan-Boltzmann constant. Finally, Qdif 
is the heat-diffusion term, defined by 



^dif 



Vth 



d 



2R dR V OR 



R- 



(10) 



and t'th(= 2W/ilI]) is the thermal diffusivity. Here, we 
use the a formula with the same a as the viscosity (see 
below). 

We employ the following viscosity parameter {a) 
prescription in a time-dependent manner (Mineshige, 
Wheeler 1989): 



a{R,t) 



ao 



H{R, t) 

'r 



(11) 



where ao is a constant, and R and H are the disk radius 
and height. From the hydrostatic relation, H = Cs/f^K, 
where Cg is the sound speed, we can rewrite equation 
(0) as 



a{R,t) 



(12) 



where 7^ is the gas constant, /i is the mean molecular 
weight, and v^{R) = R flK^R), the azimuthal velocity in 
the Keplerian disk. This implies that a is a function of 
the mid-plane or central temperature: a oc Tc{t)^^'^ . 
The boundary conditions are 

dT 

— ^ = and M = Mt (13) 
or 

at the outer edge and 
Tc = = 



(14) 



Fig. 1. Opacities from Alexander et al. (1993), Cox 
and Stewart (1970), and Cox and Tabor (1976) for 
logp = —6.0 to —8.0 at intervals of Alogp = 0.5. 



at the inner edge of the disk, where the mass-transfer 
rate from the companion star is taken to be constant in 
any particular model. Relativistic effects are not taken 
into account. 

2.3. Methods of Calculations 

The radially-dependent equations, (§)-(^, are solved 
with the equations obtained by integration of the vertical 
structure of disks. 



/-/(rc,s,i?). 



(15) 



where the function / describes the disk parameters, such 
as the fiux, Q^^^, viscous heating, Q^jg (or W), disk 
height, H, central (mid-plane) density, pc, and optical 
depth, r. 

We have taken the opacity from Alexander et al. 
(1983) for low temperatures (logTc ^ 4.0), and Cox and 
Stewart (1970) and Cox and Tabor (1976) for high tem- 
peratures (logTc ^ 4.0) with mass fractions = 0.71 
and = 0.27. The opacity for the intermediate tem- 
perature region is interpolated. The resultant opacity in 
the low-density region appropriate to our models in the 
outer disk is presented in figure 1, where the dominant ef- 
fects are designated. We plot the opacity for logpc — —6 
to —8. The disk regions in which we are interested have a 
density distribution of logpc = —6.2 to —7.7. The range, 
logTc ^ 3.25-3.55, is the region dominated by molecular 
opacity. As shown in figure 1, the effect of molecules di- 
minishes at logTc ^3.5, depending on the density. At 
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logTc ^ 3.7, H~ first becomes important. We define the 
domain between these two temperature regions as an in- 
termediate state. The partially ionized region can reach 
well above 10^ K, depending on the density. In terms of 
the effective temperature, the hot state is roughly log Tgff 

> 103-85. 

Two regimes where changes in the physical conditions 
can affect the stability of the disk have been discussed 
in previous studies. One is the partially ionized, in- 
termediate warm state centered at logTeff ^ 3.8, or, 
equivalently, logTc ^ 3.7-4.1 (Mineshige. Osaki 1985; 
Mineshige 1988; Kim et al. 1992). The other regime 
occurs in the temperature region centered at logTcff ~ 
3.4 (or logTc ~ 3.4—3.7) in a marginally optically thin 
region (Cannizzo et al. 1982; Mineshige, Osaki 1983) 
dominated by molecular opacity, such as H2O (Cannizzo, 
Wheeler 1984). These regions of variable opacity can 
lead to "metastable" or "stagnation" states in the evo- 
lution between the hot and cold stable states of the disk 
(discussed in subsection 4.2). Ionization stagnation is ex- 
plicitly displayed in models for the outburst rise by Mi- 
neshige (1988) and Kim et al. (1992). The possible sig- 
nificance of the molecular opacity to the disk-instability 
model was raised by Cannizzo and Wheeler (1984), but 
no explicit presentation has been provided. 

The basic implicit numerical code used to compute 
time-dependent disk-instability models is that given by 
Mineshige (1986; see also Mineshige, Wheeler 1989; Kim 
et al. 1992). The present disk- instability models have 
been improved compared to those presented by Mineshige 
and Wheeler (1989), in which an unreasonably small disk 
size (3.16x10^° cm) was adopted so as to avoid numerical 
instabilities and in which irradiation is ignored. These 
models represent an implicit computation of 21 radial 
zones separated by radial mesh points at fixed spacings of 
log^i? - 0.22, from logi?in = 6.5 to logi?out = 11.0. This 
zoning is adequate to reproduce the qualitative features 
of the expected outburst light curves for our adopted 
prescription for the viscosity parameter (see below). It is 
clearly not sufficient to provide a fully satisfactory model 
of a time-dependent irradiated disk. These models repre- 
sent the first study of both the radial and time-dependent 
nature of irradiated disks around black holes where the 
central star cannot provide a source of flux. They are 
presented as a step beyond the single irradiated annuli 
studied by Tuchman et al. (1990) and Mineshige et al. 
(1990b) and as a step toward the ultimate solution of 
time-dependent irradiated disks. They should not be 
viewed as full disk models, but as collections of coupled 
annuli that, nevertheless, give some qualitative insight 
into the physics of irradiated, time-dependent accretion 
disks. 



2.4- Model Parameters 

We choose a 4:Mq black hole with a companion of 
0.27Mq as being representative of A0620-00 (Marsh et 
al. 1994 and references therein). This choice, with the 
observed orbital period of 7.75 hr, gives the following 
disk model parameters: Rin = 3.6x10^ cm, binary sep- 
aration 2.24x10^^ cm, Roche lobe radius of the primary 
system 1.83xl0"'^'^cm and the disk radius lO^-'^cm (Frank 
et al. 1992; Paczynski 1977). Both the inner and outer 
boundaries of the disk are held fixed. We tested a range 
of values of -Rout, and found that the qualitative behav- 
ior is similar for outer-disk radii in the range 5xl0^''cm 
< -Rout ^ 1.3xl0^^cm, the tidal radius of the disk. We 
choose the mass-transfer rate from the companion star 
to be Mt = 3.16xlOi'5g s"^. This is larger than that 
given by McClintock and Remillard (1986), but within 
the observational uncertainty (sec the discussion in Mc- 
Clintock et al. 1995). The transferred mass from the 
companion is initially accumulated in the outermost disk 
(6x10"'^° — 10^^ cm) and then diffuses and accretes inward. 

The outburst recurrence time in these models is about 
3-4 years. This is too short for most observed systems, 
but docs not affect the qualitative thermodynamics of the 
irradiation, which is the current focus. The recurrence 
time would be longer if the transfer rate were smaller 
(see also King et al. 1997). For this relatively short 
recurrence time, the inner disk does not have time to 
return entirely to the cool state. The mass-flow rate in 
the inner disk, and hence the generated radiation, does 
drop by many orders of magnitude (about 7) from the 
peak, which is sufficient to capture the essence of the 
variable production of irradiation and how it reacts with 
the outer disk during the phases of heating and cooling. 

Given the limitation of the small number of zones, we 
have checked the basic behavior of the code by computing 
the speed of propagation of the cooling wave. Cannizzo 
et al. (1995) express the speed of the cooling front as 

V,.... = |«„(f)""5ii, (16) 

where = CgfJ is the vertical scale height at radius R, 
Gs is the sound speed, f2 is the Keplerian angular veloc- 
ity, w is the width of the cooling front and the viscosity 
parameter is given by a = aQ{H/R)^. The results of 
our numerical models are in reasonable agreement with 
this expression for our assumed a prescription and log- 
arithmic zoning. In the first 40 days after the maxi- 
mum, the velocity of the cooling wave, Vfronti is almost 
the same in both the irradiated and the non-irradiated 
models, but slightly slower in the irradiated case. At later 
times, the cooling-wave speed in the irradiated model is 
somewhat faster than that for the non-irradiated case, 
although both become very slow, less than 10"^ cm s"-'^ 
after 100 days. The overall shapes of the velocity profiles 
are rather similar in both models. 
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Although a certain range has been explored, we have 
not undertaken a systematic parameter study of these 
computationally expensive irradiated models. Varying 
Cx, ctQ, and N gives different slopes of the light curves 
in the decay and different intervals between the peaks 
of the optical activity and the inner mass flow rate. If 
the indirect irradiation efficiency parameter, Cx, is suffi- 
ciently large, Cx ^ 10~^, the models will be heated to a 
permanently hot steady state (see van Paradijs 1996). 
To approximately reproduce the observed light curve 
of A0620-00, we adopt N = 2, ^ 10^, and Cx = 
1.85x10"^. For the chosen parameters, the outbursts in 
these models always begin in the outer disk. This deter- 
mines the systematics of the overall outbiirst. In partic- 
ular, for systems with a fast rise and exponential decline, 
the optical outburst is predicted to arise well before the 
increase of the mass flow rate in the inner disk that might 
be associated with harder flux. 

We have computed several repeated outbursts in each 
model to achieve an approximate steady state. Our cur- 
rent models with irradiation show an outburst cycle of ~ 
4 yr. We did not attempt to reproduce longer recurrence 
times in this study, since our main current consideration 
concerned the effects of irradiation on the outburst, it- 
self. In the disk-instability model, the recurrence time 
scale is related to a, Mt, -Rout and the masses of the 
binary components. Mi and M2 (e.g., Cannizzo et al. 
1988; Mineshige, Wood 1989). In principle, the easiest 
way to achieve a longer recurrence time would be to de- 
crease Mt. While the value we have chosen here, Mt 
= 3.16xl0^^g s~^, is plausible, it could be a factor of 10 
smaller and still agrees with the observational constraints 
on, e.g. A0620-00. 

3. Outburst Evolution 

3.1. Light Curves 

The optical light curves are computed by assuming 
that each disk annulus radiates as a black body at tem- 
perature T(.if{R,t) and summing the emission in the V 
band over the annuli. The resultant model optical light 
curves are presented in figure 2 together with observa- 
tions of A0620-00. (The data for A0620-00 is plotted 
1.5 mag brighter than the observed value for clarity of 
presentation.) The model with both direct and indirect 
irradiation gives an optical light curve with a fast rise, 
a maximum, and a nearly exponential decline. As can 
be seen from figure 2, the direct irradiation has virtu- 
ally no effect on the light curve. The modest indirect 
irradiation we have invoked gives a slight, but noticeable 
flattening of the light curve. More severe indirect irra- 
diation would cause this flattening to be more extreme. 
For the irradiation parameter Cx ~ 10~^, the disk would 
be stabilized. The nature of this irradiation stabilization 



Fig. 2. Model light curves on the decline from the maxi- 
mum along with the optical observations of A0620- 
00 (Webbink 1978 and private communication) . The 
model optical maximum is defined as day zero. The 
data for A0620-00 is plotted 1.5 mag brighter than 
observed for clarity of presentation and their phase 
is set by assuming coincidence of the observed soft 
X-ray peak and the model optical maximum (see 
text). The model with both direct and indirect 
irradiation fits the observed optical light curve of 
A0620-00 with an assumed inclination of 60° and a 
distance of 400 pc. 



will be explored in future work. 

Note that our calculations never run into a radiation 
pressure-dominated regime nor an optically thin one, 
probably because the peak luminosity remains a slightly 
below that at which the radiation pressure dominates 
over the gas pressure. 

3.2. Mass Flow 

The model accretion rate through the inner boundary 
(Rin) of the disk is presented in figure 3. Note that since 
the mass flow rate depends sensitively on i?i as M cx Rf^f 
and since R^ is much larger than that of the marginally 
stable circular orbit (3rg) in our calculations, due to the 
crude zoning, our quiescent mass-flow rate is much larger 
than that given by the fine-mesh calculations of Ludwig 
et al. (1994). However, the qualitative behavior should 
be adequately captured by these models. 

The indirect irradiation leads to a lower Mi„ in quies- 
cence due to a higher depletion of the disk mass during 
outburst. The mass-transfer rate through the inner edge 
of the disk determines one component of the soft X-ray 
flux. In the current models, Min does not begin to rise 
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Fig. 3. Time-dependent evolution of the innermost ac- 
cretion rate (Min) for cases with no irradiation 
(short dash), direct irradiation only (long dash), in- 
direct irradiation only (dot-dash), and both direct 
and indirect irradiation (solid line). 



Fig. 4. Shadow effect for the direct irradiation. The 
hatched part in the figure signifies where a given 
zone at a given time is in the shadow of some interior 
zones that are exposed to direct irradiation. 



until 13.6 d before the primary optical peak, over two 
weeks since the start of the optical display. The peak in 
the mass flow rate occurs about 20 d after the primary 
optical maximum and about 30 d after the first rise in the 
inner mass-flow rate. Qualitatively, these models predict 
that the optical flux from the outer disk should rise be- 
fore any activity from the inner disk, soft X-rays from the 
inner optically thick disk or hard or soft X-rays from any 
coronal emission. This is consistent with the results of 
multi- wavelength observations of Nova Muscae (see Lund 
1993). A prediction of the absolute times of peak emis- 
sion in any band is more uncertain in the absence of reli- 
able models for the time-dependent power-law emission. 

3.3. The Shadow Effect 

The non-monotonic distribution of the disk height 
causes the middle portions of the disk to block the flux 
of direct radiation emitted from near the inner edge of 
the disk, and thus prevent the direct irradiation of outer 
parts of the disk. The shadowing of the outer disk is de- 
termined by the radial distribution of the opening angle, 
taiiO = H/R, assuming the radiation to arise from an in- 
ner region of small H and R. A region with a given value 
of H/R will thus be shadowed by a region at smaller R 
with larger H/R. Outer regions with small H/R in the 
wake of an inward-propagating cooling wave will clearly 
be in the shadow of inner, hot regions. 



Some systematics of the shadowing process are illus- 
trated in figure 4. The specific details reported here 
may be subject to the effects of finite zoning. The qual- 
itatively significant results arise from comparing models 
with and without irradiation. 

Figure 4 shows that the innermost disk (< 2x10^ cm) 
is always exposed to the direct irradiation. Because the 
outburst is initiated about 25 d prior to the optical peak, 
more extensive portions of the disk are exposed to di- 
rect irradiation. After maximum light, the outer disk (> 
2x10* cm) is essentially totally blocked from direct ir- 
radiation by the inner disk. With the shadowing effect, 
therefore, the existence of a mechanism which can reflect 
the radiation from the inner region is important in the 
current models if there are to be any effects of irradiation 
in the outer disk. 

3.4. Amplitude of Disk Irradiation 

The disk heating in the non-irradiated disk-instability 
model is given by the viscous heating, Q^^^ [see equation 
(^)] . In the presence of irradiation, the total heating is 
described by 

Qtot = Q+s + + ^d. (17) 

The presence of irradiation affects the disk structure, and 
hence the contribution from Q^jg, itself, is altered. The 
resultant contribution of the irradiation to the total disk 
heating is presented in figure 5 for an inner, a middle, and 
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Fig. 5. Contribution of irradiation in the disk evolu- 
tion. For illustration, we choose four zones: 5 (short 
dash), 16 (dot), 18 (dash dot), and zone 20 (solid 
line). The contribution of both direct and indi- 
rect irradiation is plotted with respect to the to- 
tal heating in the upper panel. In the lower panel, 
we present the ratio of direct to indirect irradiation. 
The direct irradiation shuts off its effect about 20 d 
after outburst maximum (day zero) due to shadow- 
ing (see figure 4). 



an outer region with the model optical maximum again 

being defined as day zero. In the inner region 10^ cm), 
the irradiation is negligible throughout the entire evolu- 
tion. The direct component is about 10 % of the indirect 
component in the inner region, but the total irradiativc 
heating does not contribute substantially compared to 
the viscous heating. In the middle of the disk (~ 10^° 
cm), the contribution of the irradiation is small, but fi- 
nite. Direct irradiation is appreciable only during the 
early decay, and is truncated by 20 d due to shadowing. 
In the outer disk (^^ 6 x 10^° cm), the contribution of 
the irradiation becomes important. Around maximum, 
the viscous heating dominates the irradiation. After the 
cooling wave propagates into the outer region, and this 
region cools, the irradiation heating, due entirely to indi- 
rect irradiation because of shadowing, dominates over the 
viscous heating. This is shown by the rapid increase in 
the ratio of irradiation heating to the total heating about 
5 d after the max;imum in figure 5. The indirect irradia- 
tion heating then decreases along with Min as the cooling 
wave propagates inward. The behavior of the outer re- 
gions under the effects of radiation is more complicated 



Fig. 6. Time-dependent evolution of the disk mid-plane 
temperature during the early decay phase (upper 
panel) and the corresponding states of the gas (lower 
panel) for the model with no irradiation. "HOT" 
denotes the hot thermal equilibrium states and 
"MOLECULE" represents the molecular stagnation 
stage (see text and figure 1). 



than reflected in figure 5. 

3.5. Global Effects of Irradiation on Decay 

The influence of mild irradiation is not restricted to a 
small modification of the light curves. Rather, it causes 
crucial changes in the local disk structure, as illustrated 
in figures 6 and 7. The upper panel of figure 6 presents 
the time-dependent evolution of the disk central temper- 
ature in the model without irradiation. This model re- 
produces the typical decay phase, as shown in numerous 
other studies of the disk instability. As the cooling wave 
propagates from the outer to inner radii, the disk makes 
the transition from the hot to cool state from larger to 
smaller radii. Note the sharp edges of the rapid transi- 
tion from the hot to cool state, which represent the cool- 
ing wave propagation. The lower panel of figure 6 shows 
the opacity domains through which the disk evolves. It 
is basically only in the hot, ionized or cool, molecular 
states with very rapid transitions between them. This 
is because in the non-irradiated model, or in the model 
with only direct irradiation and shadowing, the stagna- 
tion phase is absent. 

Figure 7 displays the time-dependent disk evolution in 
the decay phase for the irradiated disk-instability model 
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Fig. 7. Time-dependent evolution of the disk mid-plane 
temperature during the early decay phase (upper 
panel) and the corresponding states of the gas (lower 
panel) for the model with both direct and indirect 
irradiation. "HOT" denotes the hot thermal equi- 
librium states "Partial Ionization" means that the 
disk is in the metastable ionization stagnation stage, 
and "MOLECULE" (or "M") represents the molec- 
ular stagnation stage. "Intermediate State" (or "i") 
denotes the intermediate state between the molecu- 
lar and partially ionized region (see text and figure 
!)■ 



with both direct and indirect irradiation. The upper panel 
of figure 7 shows the evolution of the disk central tem- 
perature, Tc- The squares denote the time at which the 
direct irradiation is terminated at a given radius as it falls 
into the shadow of interior regions. With irradiation, the 
evolution is much more complex, as illustrated in figure 7. 
At 14.2 d, the region at ~ 4 x 10^° cm reaches the critical 
surface density for a thermal instability and begins the 
rapid cooling that characterizes the propagating cooling 
front of figure 6. The same transition is seen at ~ 2 x 10^*^ 
cm after 41.6 d and at ~ 1 x 10^" cm after 88 d. Rather 
than dropping directly to the cool state, however, these 
regions linger in an intermediate "metastable" stagnation 
state. This is obviously due to the effects of irradiation, 
since this behavior is absent in figure 6, the effects, how- 
ever, are rather subtle and involve global coupling in the 
outer regions. 

Throughout the disk evolution from quiescence to out- 
burst, the outermost disk (around 10^^ cm) stays in the 
cold, neutral state at logTes ^ 3.4, or logTc ^ 3.55. Af- 



ter maximum light, the cooling wave is initiated around 

6 X 10"'^'^ cm. The resulting outward diffusion of the sur- 
face density results in a small increase in the temperature 
at larger radii. The molecular stagnation (discussed in 
subsection 4.2) and consequent higher temperature of the 
outer regions due to the indirect irradiation affects the 
behavior of the inner regions. The regions just interior 
to the molecular-stagnation region can not simply drop 
into the cool state, because there is now some resistance 
to the outward diffusion of matter that is attendant to 
the non-irradiated cooling wave. Instead, the decline in 
temperature after the onset of the cooling instability is 
halted in the nearly constant temperature intermediate 
ionization stagnation state, as shown in figure 7. 

After the outburst maximum, the X-ray novae have 
exhibited additional features in the decay, such as "re- 
flares" 50-80 d after the maxima, the "second max- 
ima" a few hundred days later, and subsequent "mini- 
outbursts" before returning to quiescence. The reflare 
is common in many X-ray novae throughout a variety 
of wavelengths: optical (A0620-00), UV (Nova Muscae 
1991, GRO J0422-F32) and soft X-rays (A0620-00, Nova 
Muscae 1991, GRO J0422+32). As shown in figure 7, 
one portion of the cooling disk undergoes a heating in- 
stability that takes it back to the hot, ionized state, thus 
producing a transient increase in the light curve (see fig- 
ure 2). This might be associated with the reflare, since 
the optical feature is consistent with Kuulkers (1998), 
who summarized the data for the refiare in A0620-00 in 
both optical and X-rays (with some delay). It is not clear 
whether the optical feature in our models has some phys- 
ical basis or whether it is an artifact of the crude zoning. 
Our models do not provide any modulation of the inner 
mass flow and hence of the soft X-rays at the time of the 
reflare. Further studies of these issues are needed. 

4. Irradiation and Thermodynamics 

4-1- The Limit Cycle 

Figure 8 presents the actual track of the evolution 
of the surface density and mid-plane temperature T^, 
the so-called thermal limit cycle (for reviews, see Osaki 
1989; Cannizzo 1993) for various models. The top pan- 
els in flgure 8 present the evolution of the region around 
1 - 2 X 10^° cm (panel a) and around 4 — 6 x 10^° cm 
(panel c) for the case with both direct and indirect irradi- 
ation. The open circles show when the region is subject 
to direct irradiation and the closed circles show when 
the region is shadowed from direct irradiation at the in- 
dicated times. For the same two regions, the lower panels 
in figure 8 present models corresponding to no irradiation 
(long dash), direct irradiation only (short dash), and in- 
direct irradiation only (dotted line). The model with 
only indirect irradiation shows very similar evolution to 
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Fig. 8. Thermal limit cycles in the (S, Tc) plane for annuli at 2 X 10"'""cm (left panels) and 6 X lO-'-'-'cm (right panels) with both 
direct and indirect irradiation during the decay phase from the outburst maximum to the quiescent minimum state. In 
the upper (a) and (c) the open and filled circles represent epochs when the zone is and is not subject to direct irradiation, 
respectively. The outer region, at 6 X lO-'^'^cm, receives no direct irradiation due to the blocking by the inner disk and 
hence the shadowing is in effect from the outburst maximum until the disk returns to the cool quiescent state (see figure 
4). The numbered epochs correspond to times after maximum in days as follows: 1 — 0.94; 2 — 3.5; 3 - 4.7; 4 — 6.9; 5 - 14; 
6 - 15; 7 - 21; 8 - 41.6; 9 - 41.9; 10 - 44; 11 - 53; 12 - 54; 13 - 55; 14 - 57.2; 15 - 57.4; 16 - 57.5; 17 - 70; 18 - 88.2; 19 - 
88.4; 20 — 94; 21 — 130; 22 - 211. In the lower (b) and (d) the trajectory of the annulus evolution is illustrated for three 
cases: non-irradiated (dotted line), directly irradiated only (dashed), and indirectly irradiated only (solid line) models. 



the full model. There are three distinctive stages to the 
full model: 

(1) In the region from 1 — 2 x 10^" cm, the temperature 
stops declining and then climbs back up to the hot state 
again (figure 8a, time steps 9—14). 

(2) On the second decline from the hot state the disk tem- 
perature stagnates at logTc ~ 3.8—4.1 (figure 8a, time 
steps 16—19) in the partially ionized region (see figure 7 
and subsection 4.2), prior to the final downward transi- 
tion to the cool disk. We call this phase an "ionization 
stagnation." 



(3) In the outer region, ^ 4 — 6 x 10^" cm, the disk shows 
a similar stagnation, but it occurs at logTc ^ 3.4—3.6 
(figure 8c, time steps 3—7) in the molecular opacity- 
dominated region (see figure 7 and section 4.2). We call 
this phase "molecular stagnation." 

Figure 8 shows that the stagnation behavior on the de- 
cline is not present in the non-irradiated or only directly 
irradiated models. Even with modest indirect irradia- 
tion this stagnation behavior can affect the outburst and 
cooling of the disk. 
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4-2. Stagnation Phenomena 

The disk is in thermal equilibrium when the viscous 
heating, Q^^^ [equation (H)], balances the radiative cool- 
ing, Q-^^ [equation (|)], or = Q'^^. In the thin disk 
approximation (e.g., Shakura, Sunyaev 1973), the emer- 
gent heat flux radiated vertically in the plane-parallel 
approximation is given by F(z) = 4(7Tc^/3t. Since r — 
Ek for the Rosseland mean opacity k and the surface 
density is S = 2pH for density p and disk height H, one 
can write 



ST," 



(18) 



Since the sound speed is Cg = {TZT,/ pY^"^ , where TZ is 
the gas constant and /i is the molecular weight, and H 
= Cs/r2, the energy equation in thermal equilibrium be- 
comes 



3 ^(nr, 

2 V 



3 P iacT," 
it = —ail — = 

2 p 3kS] 



(19) 



where a = ac/4 and k = k{Tc, Pc)- Note that, in the 
disk-instability model which we present here, a is also 
a function of [see equation (|ll])]. In various physical 
regimes, the opacity can be expressed as 



K = KqPc 



(20) 



where kq is a constant and to > 0, but n can be negative 
or positive (see figure 1). 

The possibility of stagnation in a partially ionized re- 
gion was first pointed out by Meyer-Hofmeister (1987) 
and elaborated by Mineshige (1988). For the case of a 
departure from thermal equilibrium, the energy equation, 
equation (Bh, can be written as 



2 dt 



(21) 



where Cp is the specific heat at constant pressure, and 
Q"*" and Q~ are the total heating and cooling rates. The 
rate of change of the temperature can be decreased by in- 
creasing the specific heat, or by decreasing the departure 
from thermal equilibrium, Q"*" ~ Q~, or a combination of 
both. The essence of the higher temperature stagnation 
is that the specific heat increases in the regions of partial 
ionization and if the structure is not too far from thermal 
equilibrium, the thermal time scale is increased, leading 
to a "metastable" intermediate temperature state. In the 
rise, the partially ionized region is near the "knee" in the 
instability curve, where ^ and the effect can be 
pronounced. 

Mineshige (1988) showed that the onset of the heat- 
ing can consist of a "warming wave" which raises the 
disk from the cold state to this intermediate metastable 
state. This state is ultimately unstable and a full heat- 
ing wave causing the transition to the hot state finally 



ensues. In the decline, the partial ionization state is fur- 
ther from the equilibrium curve, and although the spe- 
cific heat increases, the cooling term, Q^, dominates the 
heating term, and there is little effect of the specific heat 
as the disk drops into the cool state. As a result, stagna- 
tion has not been manifested in the decline in previous 
disk-instability models. As the current models illustrate, 
irradiation can enhance the heating term, Q+, in the de- 
cline, thus restoring an approximate balance, Q"*" ~ Q~. 
In this situation, the effect of an increase of Cp can be 
seen as "stagnation" in the decline as well. This can 
cause the disk to linger in the intermediate-temperature 
state. It eventually makes a downward transition into 
the cool state. 

At logTc ^ 3.5—3.6, the molecular absorbers become 
important and at logTc ^ 3.4, H2O dominates. The low- 
temperature opacity is dominated by grains as grain con- 
densation begins at logTc ^ 3.25. These various effects 
change the slope of the opacity at about logTc — 3.25, 3.4 
and 3.5—3.6 (figure 1). As a result, when the disk passes 
through logTc = 3.6, the opacity suddenly increases with 
decreasing temperature [see equation (p^)]. As a result, 
the thermal time scale for the downward transition will 
be increased. 

5. Discussion and Conclusions 

The purpose of this paper was to investigate some 
of the basic time-dependent effects of irradiation. We 
show that even relatively mild irradiation from the inner 
disk can affect the overall disk evolution. Mild irradia- 
tion can affect the limit cycle by causing portions of the 
disk to linger in the intermediate-temperature metastable 
"stagnation state." In contrast to the non-irradiated disk- 
instability models, the irradiated models show a wave of 
cooling to the stagnation state that propagates inward in 
advance of cooling to the low, molecular-dominated state. 
It is not clear that this is simply a result of heating of 
the disk, since the non-linear interplay of the zones seems 
to have an important role. The issue of how a particu- 
lar portion of the disk responds to irradiation by tracking 
along the metastable state is a rather subtle thing involv- 
ing complex interactions among the different portions of 
the disk. Due to shadowing, the effect of direct irradi- 
ation terminates as each annulus of the disk enters the 
stagnation state. The net result of the irradiation is a 
slower decay of the light curve. 

When an outburst is initiated, the outer disk is exposed 
to direct irradiation, but that radiation is feeble until the 
heating wave reaches the inner disk. Soon thereafter the 
outer disk is shadowed by a swelling of the middle of the 
disk. Direct irradiation does not play an important role 
in the behavior of the outer disk either during and after 
the maximum of the outburst, but can affect the middle 
of the disk during the outburst. Since it is by assumption 



No. 3] 

not influenced by shadowing, the indirect irradiation af- 
fects the overall disk evolution, including the outer disk, 
throughout its evolution. 

As stated in the introduction, this is the first study of 
time-dependent irradiation of black-hole accretion disks. 
For this reason and technical issues related to code insta- 
bility, we have presented models with rather crude zon- 
ing. These models are a major step beyond the one-zone 
models of Tuchman et al.(1990), but only a step toward 
the goal of self-consisent, resolved irradiated models. Be- 
side the rather crude zoning, other restrictions are rather 
simplified models for the irradiation that make no at- 
tempt to do realistic radiative transfer and, of course, 
continued uncertainty about the physics of disk viscos- 
ity and the role of disk evaporation. In the context of 
these general limitations and uncertainties, these models 
do give some insight into the issues of time-dependent ir- 
radiation. They incorporate sufficient physics to capture 
the non-monotic nature of the disk profile, a critical issue 
in considering shadowing. Given this context, there are 
some qualitative aspects of these models that are reason- 
ably trustworthy: 

• The direct irradiation of the disk by inner layers 

where the soft X-ray flux is generated is found to 
have only a small effect on the outer disk because of 
shadowing. 

• Mild indirect irradiation that flattens, but otherwise 

docs not affect the light curve substantially, still has 
interesting non-linear effects on the structure of the 
disk as heating and cooling waves propagate. 

Any quantitative results (e.g., front propagation velocity, 
the detailed shape of the light curves, etc) may change 
when fine-mesh calculations are performed. Specifically, 
the calculated reflare-like optical events may or may not 
be a calculation artifact. We need fine-mesh calculations 
to confirm the reflare-like behavior of the irradiated disk, 
but the basic phenomenon of radiation-induced stagna- 
tion allowing a temporary return to the hot state is phys- 
ically plausible. It is not clear how such a phenomenon 
could lead to a modulation of the soft X-rays, a defining 
feature of the observed rcflares. 

Saito (1989) considered the case of a central star that 
illuminated the disk using a simple one-zone model for 
the vertical disk structure, and thus did not properly 
consider the non-(thcrmal)-equilibrium states of the disk. 
The irradiation from a central source is larger than what 
we consider hero by a factor of H/ R (cf. subsection 2.1, 
van Paradijs 1996). In addition, Saito neglected a term 
equivalent to [dhiH/d\nR- 1] = {R/ H)[d{H / Rf / dR] 
[see equation (3)], that measures the shape of the disk 
profile. He took the angle of incidence of the radia- 
tion, ~ H/R, and the fiux generated by accretion to be 
constant. Despite these differences and approximations, 
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Saito's results are qualitatively similar to ours. With a 
large accretion rate (10^^ g s^^) and the various approx- 
imations that enhance the irradiation, Saito found that 
the inner disk was always maintained in a permanently 
hot state. The outer disk, however, underwent a thermal 
instability that began in the outer part of the disk. Saito 
notes that prior to the outburst the outer disk is in the 
shadow of the inner disk. 

The disk instability due to the ionization of hydrogen 
and helium remains the most plausible cause of the out- 
burst of the black-hole candidate X-ray novae. For the 
orbital periods and mass-transfer rates inferred from ob- 
servations, the disk is predicted to be unstable. A steady 
state is very unlikely. Because the direct component of 
irradiation comes from a surface of the disk that is nearly 
in the plane of the outer disk, and hence presents a small 
solid angle, disks around black holes receive relatively lit- 
tle direct irradiation (Cannizzo 1994; King et al. 1997; 
Cannizzo 1998; Dubus et al. 1999). Gradients in opacity, 
the viscosity, the effect of the inner boundary condition, 
and of time-dependent structure, all serve to diminish the 
direct irradiation even further. Whether a source of "in- 
direct" irradiation as we have modeled it here can have 
a more severe effect, including stabilizing the disk (van 
Paradijs 1996), requires further investigation. 

The value of the indirect irradiation parameter, Cx 
10~^, was chosen in conjunction with other parameters in 
these models to reproduce the approximate slope of the 
optical light curve of A0620-00. We find that this level of 
irradiation does not significantly affect the luminosity at 
maximum, but does make the disk brighter by about 0.5 
mag at about 50 d after maximum. Cannizzo (1998) has 
estimated that A0620-00 might be brighter by about 1 
mag at maximum light than his standard non-irradiated 
models. This difference may be within the model un- 
certainties, but, if ascribed to irradiation, it could be 
accounted for with a modest increase in Cx- 

We have shown that even modest levels of irradiation 
can influence the thermodynamics of the disk and affect 
its evolution. It is very important to consider indirect as 
well as direct irradiation of the disk, since the effect of 
the latter is severely constrained by shadowing. 

Finally, we can consider how our results might change 
if the quiescent disk has a hole in the middle because 
of gas evaporation and resultant advection-dominated 
flow. Such a hole will change, for example, the quiescent 
surface-density distribution, thereby affecting the igni- 
tion time and radius (Hameury et al. 1997). Also, the 
disk should expand during quiescence because of continu- 
ous input of angular momentum by the incoming stream, 
while the mass increase will be modest (Mineshige et al. 
1998). Once an outburst occurs, however, the entire disk 
evolution will be totally controlled by the propagation of 
the heating front, except at the very early rise phase, so 
that our conclusions will not change signiflcantly. 
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